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Comment on: Reply to comment on 'Perfect imaging without negative refraction' 

Paul Kinsler* and Alberto Favaro 
Blackett Laboratory, Imperial College, Prince Consort Road, London SW7 2AZ, United Kingdom. 

(Dated: Friday 28 th October, 2011) 

Whether or not perfect imaging is obtained in the mirrored version of Maxwell's fisheye lens is 
debated in the comment/reply sequence [1, 2] discussing Leonhardt's original paper [3]. Here we 
show that causal solutions can be obtained without the need for an "active localized drain" , contrary 
to the claims in [2]. 
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R.J. Blaikie (RJB) notes in [1] that that Ulf Leon- 
hardt's (UL) setup in [3] incorporates an "active localized 
drain" at the image point. It is this drain, modeled as a 
phase-delayed mirror image of the source, which provides 
the sub-wavelength detail of the source's image. RJB 
showed that a steady-state numerical simulation of the 
fields without the drain did not show the perfect super- 
resolution image of the source. In response, UL noted 
[2] that steady-state solutions neglect causality and that 
inclusion of a sink solves the problem of energy buildup. 

We now address both of UL's concerns about RJB's 
simulations using a time-domain numerical solution with 
a source only active for a finite time. In figs. 1 & 2 
the results from such a simulation, calculated using the 
open source MEEP [4] implementation of FDTD [5] , are 
shown. In FDTD, the real valued electric and magnetic 
fields are defined over space, and an algorithm is applied 
which propagates these fields, step by step, forwards in 
time. It is thus explicitly causal, and entirely indepen- 
dent of any decomposition into plane waves or modes. 



therefore correctly tests for the geometric response and 
achievable spatial resolution, irrespective of the fact that 




FIG. 1: Snapshot of the electric field E z from a simulation of 
the cylindrical mirrored fisheye, as the light first reaches the 
image point. Parameters are taken from RJB but with zero 
losses: n(r) = 2/(1 + (r/V ) 2 ) for r = Wfim, / = lOOTHz 
(i.e. Ao = Sum). 

The point-like source used is independent of all other 
properties of the simulation[12]. We follow UL and RJB 
and use a frequency independent refractive index pro- 
file, so there is no temporal dispersion. Our simulation 
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FIG. 2: Time history of the simulation in fig. 1 showing the 
field E z along the x-axis. The intensity FWHM of the source 
and image(s) are given, and can be compared to either Ao or 
the local wavelength (Ao/n ~ 2.3/xm) at the image point. Fur- 
ther bounces gradually degrade the image quality, and more 
so for briefer source durations. 



that the image is not as sharp as the source, and matches 
the steady-state results of RJB [1]. 

We do not include an active drain in our simulations 
because we believe it unlikely to form a part of any ac- 
tual device; e.g. if this mirrored fisheye replaced the 
elliptical cavity used in lamp-pumped lasers. Active 
drains often appear in the literature when systems with 
sources are designed using folded-space or mirror-imaged 
transformations: e.g. the transformation optics slab-lens 
[6]. However some authors insist that such active drains 
are unphysical and/or mathematically ambiguous [7-10]. 
RJB was able to replace the drain in his steady-state 
simulation with a carefully-phased source [1], but this 
will fail in general for both time domain simulations and 
physical devices [11]. Whatever the method, we consider 
achieving super-resolution by such means to be of little 
utility, since it requires a customised element to enhance 
and 'image' the field at each precisely tuned pixel. 

In summary, despite the claims in [2], causality does 
not require the presence of an active drain, irrespective of 
whether or not an active drain might be otherwise useful. 
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Appendix: Sample MEEP control file 

Maxwell fisheye in 2d 

Paul Kinsler / Imperial College London, May 2010 



Treat sizes as if in microns (hence radius 10 ==> radius lOum) 



; Fisheye sizes and parameters 

(def ine-param n 2) ; base index of fisheye 
(def ine-param w 1) ; width of waveguide 

(def ine-param r 10) ; inner radius of ring 



; Computational resolutions etc 

(set-param! resolution 20) ; cells per unit size 

; NB: 10 is ok for a couple of bounces 
; but increase it for longer sims 

(def ine-param pad 2) ; padding between waveguide and edge of PML 

(def ine-param dpml 2) ; thickness of PML 

; Here we make sxy non-integer so that the simulation will have an 

; odd-number x and y cell count; this is because I want an exact on-axis 

; readout of fields etc without having to average array elements. 

> 

(define sxy (+ (/ 1 resolution) (* 2 (+ r w pad dpml)))) 
(set! geometry-lattice (make lattice (size sxy sxy no-size))) 



; Refractive index profile functions 

; to calculate the local refractive index as it varies with position. 

3 

; DEFINE frr = sqrt(x~2 + y~2) 
j 

(define (frr p) 

(sqrt (+ (* (vector3-y p) (vector3-y p)) 
(* (vector3-x p) (vector3-x p)) 

) ) ) 

; DEFINE fdivisor = 1 + rr~2/r~2 
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(define (fdivisor rr) 

(+ 1 ( / (* rr rr) (* r r) ) 

) ) 

; DEFINE refindex = n / fdivisor = n / ( + rr~2/r~2) 

(define (refindex rr) 
(/ n (fdivisor rr)) 

) 



; DEFINE eps = refindex~2 = the relative permittivity 

(define (eps rr) 

(* (refindex rr) (refindex rr)) 

) 

; DEFINE f medium (x,y) 

(define (fmedium p) 
(make medium 

(epsilon (eps (frr p))) 

) ) 



; Create the Fisheye structure 

; Create a metal ring to form the mirror using two overlapping cylinders 
; - later objects take precedence over earlier objects, so we put the 
; outer cylinder first, and the inner (air) cylinder second. 

; The first three blocks are just to add some guides for the 

; eye locating the source and image points — 

; they have no effect inside the mirrored fisheye 

(set! geometry 
(list 

(make cylinder (center 0) (height infinity) 
(radius (+ r w) ) 

(material perfect-electric-conductor) 
(material perfect-magnetic-conductor) ) 
(make cylinder (center 0) (height infinity) 
(radius r) 

(material (make material-function 

(material-f unc fmedium) ) ) ) 

; Some objects we might like to try putting at the image point: 

; (make cylinder (center -7.20 0) (height infinity) 
; (radius 0.02) 

; (material metal) 

; ; (material perfect-magnetic-conductor) 

; ; (material (make dielectric (epsilon 1.000001) 
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(E-polarizations 
(make polarizability 

(omega 0.01) (gamma 1) (sigma 1.0e5))))) 



) ) 



; Set up the PML at the simulation boundaries 

(set! pml-layers (list (make pml (thickness dpml)))) 



SOURCES: Put a single point source on the y-axis at x=7.20um 



(def ine-param ffq (/ 13)) 

(def ine-param f delta (* ffq 20.25)) 
(def ine-param f duration (/ 6 ffq)) 



scaled inverse wavelength of 3um 
a time delay that might be useful 
a on/off time scale for switching 



(set! sources 
(list 

Finite-time source at x=7.20um with gradual switch on/off 
to avoid generating high frequency components that will 
clutter the output (but would propagate just the same 
in this dispersionless system) . 

(make source 

(src (make continuous-src (frequency ffq) (width (* 2 ffq)) 

(end-time fduration))) 
(component Ez) (center 7.20 0.00) 

) 

Trial (pseudo-drain) source 
(make source 

(src (make continuous-src (frequency ffq) (width (* 2 ffq)) 

(amplitude -1) 
(start-time fdelta) 
(end-time (+ fduration fdelta)))) 
(component Ez) (center -7.20 0.00) 

) 



Run the simulation 



(run-until 132 

(at -beginning output-epsilon) 
(to-appended "I" 
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(at-every 0.125 (synchronized-magnetic output -t ot -pwr ) ) 

) 

(to-appended "ez" 

(at-every 0.125 output-ef ield-z) ) 

) 
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